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Abstract 

It  has  been  frequently  observed  that  the  backward  differentiation 
approximation  to  the  solutions  of  £x'  +  Fx  =  f  can  fail  to  converge 
even  pointwise  in  an  initial  boundary  layer.  This  note  shows  that  the 
approximations  converge  in  a  distributional  sense  even  if  the  exact 
solution  is  also  distributional. 


1  Introduction 

In  a  fund2unental  series  of  papers  [3],  [4],  [5],  Cobb  investigated  the  distri¬ 
butional  solutions  of  the  linear  time  invariant  descriptor  system 

Ex  +Fx  =  f  (1) 
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and  showed  that  the  distributional  solutions  were  the  limits,  in  a  distribu¬ 
tional  sense,  of  the  solutions  of  the  regularized  system  {E  +  €F)xf  -I-  Fz  =  / 
as  e  — »  0.  These  impulsive  solutions  were  also  studied  via  transfer  functions 
by  Verghese  [9].  One  of  the  first  class  of  numerical  methods,  other  then 
the  reduction  of  the  pencil  \E  -f-  F,  to  be  applied  to  (1)  were  backward 
differentiation  formulas  (BOF)  [7],  [8].  In  [8]  it  was  shown  that  the  BDF 
approximations  converge  to  the  true  solution  outside  of  an  initial  boundary 
layer  of  nonconvergence. 

The  simplest  of  the  BDF  methods  is  the  implicit  Euler’s  method  with 
fixed  stepsize  h  which  is 

=  /,+,  (2) 

where  =  fo  +  lA,  /  =  f{ti),  and  Xj  is  the  approximation  of  x(<,).  In  [2]  it 
is  observed  that  the  boundary  layer  errors  look  a  lot  like  approximations  of 
distributions.  This  note  will  show  that  the  BDF  approximations  actually 
converge  in  a  distributional  sense  on  the  entire  interval  of  definition  to  the 
solution  of  (1)  even  if  th^t  solution  is  distributional. 

2  Results 

We  assume  that  E,  F  are  n  x  n  matrices  and  that  XE  -I-  F  is  a  regular 
pencil  so  that  det(AF  +  F)  ^  0.  We  abo  assume  that  E  is  singular,  and,  in 
order  to  avoid  trivial  special  cases,  that  the  index  of(l)islr>2.  That  is, 
the  order  of  the  pole  at  infinity  is  at  least  two.  Equivsdently,  the  solutions 
of  (1)  will  involve  Jt  —  1  derivatives  of  the  forcing  function  /.  We  assume 
that  /  is  at  least  k  times  continuously  differentiable.  In  order  to  establish 
convergence  for  general  fixed  step  BDF  methods,  it  sufBces  to  consider  (2). 
Notation  is  simplified  by  taking  the  initial  time  to  to  be  zero. 

Given  a  sequence  {x,}~o  identify  it  with  the  function 

f;  -  i/i)  -  ff(<  -  (i  +  i)A)] 

1=0 

where  H  is  the  unit  step  (Heaviside)  function.  If  {x,}  is  given  by  (2),  then 
it  gives  an  0(A)  approximation  to  the  true  solution  x(f)  on  the  interval 
[FA  T]  for  fixed  T  [8].  We  are  interested  in  what  happens  on  [0  AA]. 
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By  using  the  Kronecker  structure  of  the  pencil,  and  the  well  known 
theory  of  singular  systeirvs  [1],  [2],  it  suffices  to  consider  the  special  case 


Nx*  +  x  =  f 

where  N  is  k  x  k  nilpotent  Jordan  block 

0  1  •  0 
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Then  (2)  becomes 

{N  +  hl)xi+i  =  Nxi  +  hfi+i 

We  first  consider  the  associated  homogeneous  equation  for  (3) 


(4) 


Nx'  +  1  =  0,  z(0)  =  lo  (5) 

ifhe  solution  of  (5)  is 

k~2 

x  =  Y.i-iy6^'^N^^'^Xo  (6) 

t=0 

where  is  the  itb  derivative  of  the  delta  function  (5(f).  The  backward 
Euler  approximation  for  (5)  given  by  (4)  is 


x.+i  =  {N  +  hiyK^x,  =  [{N  +  /i/)-^V]‘+'io 


which  is  zero  for  i  +  1  >  Let  p  =  l/h.  Then 

■  0  p  -p2  •  (-i)V''-'  ■ 

0  0  p 

(N  +  hi)-' y  =  \  .  .  _p2 

p 

0  ■  •  0 


(7) 


and  for  1  <  i  <  Jb  —  1, 


P‘ 


[(^’  +  hI)-^N]'  = 


where  the  ft + j  +  l)th  superdiagonal  of  (8),  for  0  <  j  <  k—l  —  i  has  entries 


Using  (8)  and  comparing  (7)  to  (6)  and  equating  the  coefficients  of  like 
coordinates  of  Xq,  we  find  that  we  need  to  show  that 
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ss  h  —*0*  where 


4^'  =  ( j- )  m  -  a + i)A)  -  m  -  u + 2)A)i 

and  (9)  is  taken  in  the  weak  distributional  sense.  That  is, 

Urn  I"  z^\t)git)dt  =  6^\t)g{t)dt  =  (-1)*^<’)(0)  (10) 

for  every  infinitely  differentiable  test  function  g  on  [0  oo).  Let  w'{t)  -  g(t) 
so  that 


=  (  ]  j[MU  +  2)b)-wiU  +  l)h)] 
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Define  the  operators  5,  A  by  Sxti  =  tx,+i  and  Arz,  =  rz,+i  —  iz,  so  that 
A  s=  5  —  Then  (11)  becomes 

=  i 


-A(/  -  5)'5tz; 


(-ly 


But  converges  to  tz;^*"'"*)  s=  [6]  amd  thus  (9)  holds. 

Now  consider  the  nonhomogeneous  problem  (3).  By  linearity  the  so¬ 
lution  of  (3)  is  made  up  of  a  smooth  solution  on  [0  T]  and  (possibly)  an 
initial  impulse  satisfying  (3)  with  /  =  0.  The  preceding  argument  shows 
that  the  approximation  for  the  distributional  part  converges.  It  sufiices 
then  to  consider  the  approximation  of  the  smooth  solution.  Let  xqt  the 
exact  initial  value  of  the  smooth  solution  of  (3).  Simple  examples  show 
that  for  the  first  k  steps  the  error  in  using  (4)  can  be  unbounded  as  h  goes 
to  0  even  if  Xq  =  Xqt-  We  wish  to  show  that,  in  fact,  the  error  goes  to 
zero  weakly  in  the  distributional  sense.  By  iterating  (4)  backwards,  it  can 
be  shown  that  there  is  am  initial  condition  xoh  such  that  taking  xq  =  xoh 
leads  to  a  solution  of  (4)  which  gives  a  uniformly  0{h)  approximation  to 
the  smooth  solution  of  (3)  on  [0  T]  and  also 

Xqk  =  XoT  +  h(p{h)  (12) 

where  ^  has  a  series  expansion  <t>(h)  =  -f  Now  let  xf'  be 

the  solution  of  (4)  with  xq  =  xo/, ,  while  xj^'  is  the  solution  of  the  aissociated 
homogeneous  equation  (7)  with  xq  =  xor  and  x^*"'  is  the  solution  of  the 
associated  homogeneous  equation  with  xq  =  <p,-  Then 

x!°’  =  xl^  +  /i£;/i'xp”’-HC>(/i) 
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But  is  a  uniformly  0(h)  approximation  of  the  unique  smooth  solution 
of  (3).  From  the  argument  used  to  prove  (9)  we  have  that  the  converge 
to  distributions,  so  that  converges  to  zero  as  a  distribution  as  h 

goes  to  zero.  In  summary,  we  have  shown: 

Theorem  1  Suppose  ihai  \E  -\-F  is  a  regular  pencil  and  Xo  is  an  arhiiray 
initial  condition.  Let  {xj}  be  the  backward  Euler  approximation  using  (2). 
Then  this  approximation  converges  weakly  in  the  distributional  sense  to  the 
distributional  solution  of  Ex'  +  Fx  =  /,  x(io)  =  xq. 

3  Comments 

We  consider  the  main  result  of  this  paper  to  be  primarily  of  theoretical 
interest.  However,  it  is  interesting  to  note  that  if  one  has  the  system  (1) 
and  is  interested  in  the  possible  impulsive  behavior,  then  impulsive  behavior 
can  be  modeled  by  using  (2).  This  simulation  is  much  quicker,  and  easier 
to  program,  then  a  code  to  compute  the  pencil  decomposition.  Abo  note 
that  if  the  quantities  of  interest  involve  weighted  integrab  of  the  solution, 
then  these  quantities  can  abo  be  estimated  using  backward  differentiation 
formulas. 
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